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Q I Abstract 

' We show that minimal SO(IO) Grand Unification models where 

the fermions have Yukawa couplings to only one (complex) 10 and one 
^ , 126 of Higgs scalars lead to a very predictive neutrino spectrum. This 

H ' comes about since the standard model doublet contained in the 126 of 

Higgs (needed for the see-saw mechanism) receives an induced vacuum 
expectation value at tree-level, which, in addition to correcting the 
bad asymptotic mass relations = mg and = m^, also relates the 
Major ana neutrino mass matrix to observables in the charged fermion 
sectors. We find that (i) the i^e — mixing angle relevant for the solar 
neutrinos can be considerably smaller than the Cabibbo angle and lies 
in the range sin^e/x = — 0.3, (ii) ~ Vt mixing is singer — 3|T4rf| ~ 
0.05, (iii) the v^^—v^ mixing angle is large, sin^^,- ~ 3| Vcbl = 0.12—0.16, 
and (iv) vap^jmy^ > 10^, implying that — oscillations should be 
accessible to forthcoming experiments. 
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It is quite possible that the deficit of solar neutrinos reported in the 
Chlorine/ Kamiokande,^ SAGE,"^ and GALLEX"^ experiments is an indica- 
tion that the neutrinos have masses and mixings very much like the quarks. 
The observed deficit can be explained in terms of neutrino oscillations in two 
different ways: (i) long wave length vacuum oscillation,^ and (ii) resonant 
matter oscillation (the Mikheyev- Smirnov-Wolfenstein (MSW) effect®). As- 
suming a two-ffavor V(. — v^ oscillation, in the former case, the neutrino masses 
and mixing angle should satisfy^ (at 90% CL) Am^ ~ (0.5 to 1.1) x 10"-*^° eV'^ 
and sin^2^e/^ ~ (0.75 to 1). In case of MSW, on the other hand, there are two 
allowed windows that fit all of the experimental data^: (a) the small mixing 
angle non-adiabatic solution, which requires Am^ ~ (0.3 to 1.2) x 10^^ eV'^ 
and sin^2^^e^ ~ (0.4 to 1.5) x 10^^, and (b) the large angle solution with 
Am^ ~ (0.3 to 5) X 10~^ eV'^ and sin226'e^ ~ (0.5 to 0.9). In all these cases, 
barring an unlikely scenario of near mass degeneracy among neutrinos, either 
Vy^ or Vr should have mass in the (10^^ to 10^^) eV range. 

A well-known and elegant explanation for the origin of such tiny neutrino 
masses is the see-saw mechanism,^ wherein the light neutrino masses scale 
inversely with the B — L breaking Majorana mass M of the right-handed 
neutrino: rrii, ~ m'^/M, m being the neutrino Dirac mass. The smallness 
of rriiy is then understood in terms of the heaviness of the B — L breaking 
scale. The solar neutrino puzzle indicates that the B — L scale is in the 
(10^2 _ xQie) Gel/ range. 

All of the observations above, viz., non-zero neutrino masses, the see-saw 
mechanism, and a high B — L scale, fit rather naturally in grand unified mod- 
els based on the gauge group 5*0(10)^. In its non-sup ersymmetric version, 
experimental constraints from proton life-time and the weak mixing angle 
sin^^VK require that 5*0(10) breaks not directly into the standard model, but 
at least in two steps. In a two-step breaking scheme, the left-right symmet- 
ric intermediate scale is around 10^^ GeV.^'' In supersymmetric 5*0(10), on 
the other hand, there is no need for an intermediate scale, 5*0(10) can break 
directly to the standard model at around 10^® GeV. 

To confront 5*0(10) models with the solar neutrino data, one must make 
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precise predictions of the neutrino masses and mixing angles. This requires, 
however, detailed information of the Dirac neutrino mass matrix as well as 
the Majorana matrix. In grand unified theories (Guts), it is possible to relate 
the quark masses with the lepton masses. For example, in SU{5) Gut, the 
charge —1/3 quark mass matrix is equal to the charged lepton matrix at the 
unification scale. Most SO (10) models retain this feature; in addition, they 
also relate the neutrino Dirac mass matrix to the charge 2/3 quark matrix. 
However, there is no simple way, in general, to relate the heavy Majorana 
matrix M to the charged fermion observables. 

The purpose of this Letter is to show that in a class of minimal 5*0(10) 
models, in fact, not only the Dirac neutrino matrix, but the Majorana ma- 
trix also gets related to observables in the charged fermion sector. This leads 
to a very predictive neutrino spectrum, which we analyze. We use a sim- 
ple Higgs system with one (complex) 10 and one 126 that have Yukawa 
couplings to fermions. The 10 is needed for quark and lepton masses, the 
126 is needed for the see-saw mechanism. Crucial to the predictivity of 
the neutrino spectrum is the observation that the standard model doublet 
contained in the 126 receives an induced vacuum expectation value (vev) 
at tree-level. In its absence, one would have the asymptotic mass relations 
rrib = rrir, nis = m^, rrid = rrie, as in minimal SU{5). While the first relation 
would lead to a successful prediction of at low energies, the last two are 
in disagreement with observations. The induced vev of the standard doublet 
of 126 corrects these bad relations and at the same time also relates the Ma- 
jorana neutrino mass matrix to observables in the charged fermion sector, 
leading to a predictive neutrino spectrum. 

We shall consider non-Susy 5*0(10) breaking to the standard model via 
the SU{2)l X SU{2)r x SU{A)c = G224 chain as well as Susy-5O(10) breaking 
directly to the standard model. Our predictions on the neutrino mass ratios 
and the mixing angles are essentially independent of the chain of descent, it 
only affects the overall scale of neutrino masses. 

The breaking of 5*0(10) via G224 is achieved by either a 54 or a 210 
of Higgs. The 210 also breaks the discrete D-parity,^^ the 54 preserves it. 
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D-parity is a local discrete Z2-subgroup of 5*0(10), under a fermion field 
/ transforms into its charge conjugate f^. Breaking of D-parity at the Gut 
scale makes the see-saw mechanism natural. The second stage of symmetry 
breaking goes via the 126. Finally, the electro-weak symmetry breaking 
proceeds via the 10. Note that in Susy-5'O(10), the first two symmetry 
breaking scales coalesce into one. 

Let us turn attention to the fermion-Higgs Yukawa couplings of the 
model. Denoting the three families of fermions belonging to 16-dimensional 
spinor representation of 5*0(10) by ipa-, a = 1 — 3, the complex 10-plet of 
Higgs by if, and the 126-plet of Higgs by A, the Yukawa couplings can be 
written down as 

Ly = hab^a^bH + fab^aA'^ + H.C. (l) 

Note that since the 10-plet is complex, one other coupling ipa'ipbH is allowed 
in general. In Susy-5'O(10), the requirement of supersymmetry prevents such 
a term. In the non-Susy case, we forbid this term by imposing a U{1)pq 
symmetry, which may anyway be needed in order to solve the strong CP 
problem. 

The 10 and 126 of Higgs have the following decomposition under O224: 

126 ^ (1,1, 6) + (1,3, 10) + (3,1, To) + (2, 2, 15) 
10 ^ (1,1, 6) + (2, 2,1) (2) 

Denote the (1, 3, 10) and (2, 2, 15) components of A(126) by and E re- 
spectively and the (2, 2, 1) component of H{10) by The vev < A^j >= 
vr ~ 10^^ GeV breaks the intermediate symmetry down to the standard 
model and generates Majorana neutrino masses given by fvR. $ contains 
two standard model doublets which acquire non-zero vev's denoted by Ku 
and Kd with Ku,d ~ 10^ GeV. Ku generates charge 2/3 quark as well as 
Dirac neutrino masses, while Kd gives rise to —1/3 quark and charged lepton 
masses. 

Within this minimal picture, if Ku, Kd and vr are the only vev's contribut- 
ing to fermion masses, in addition to the SU{5) relations = m,-, = 
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'm-fi, iTT-d = iTT-e, eq. (1) will also lead to the unacceptable relations rriu '■ Trie '■ 
rrit = rud '■ ttLs : rrih. Moreover, the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix will be identity. Fortunately, within this minimal scheme, we 
have found that there are new contributions to the fermion mass matrices 
which are just of the right order of magnitude to correct these bad relations. 
To see this, note that the scalar potential contains, among other terms, a 
crucial term 

Vi = XAAAH + H.C. (3) 

Such a term is invariant under the U{1)pq symmetry. It will be present in 
the Susy SO{10) as well, arising from the 210 F-term. This term induces 
vev's for the standard doublets contained in the E multiplet of 126. The vev 
arises through a term AjiAjiE^ contained in Vi. This induced vev can also 
be seen by analyzing the one-loop graph involving neutrinos which generates 
a divergent contribution to such a term. 

We can estimate the magnitudes of the induced vev's of E (denoted by Vu 
and Vd along the up and down directions) assuming the survival hypothesis 
to hold: / 2 \ 

Suppose Mu ~ 10^^ GeV, Mi 3 x 10^^ GeV and Ms ~ 10^^ GeV, consis- 
tent with survival hypothesis, then f„ and Vd are of order 100 MeV, in the 
right range for correcting the bad mass relations. We emphasize that there 
is no need for a second fine-tuning to generate such induced vev's. In the 
Susy version, since there is no intermediate scale at all, the factor (vji/M^) 
is not a suppression, so the induced vev's can be as large as k^j- 

We are now in a position to write down the quark and lepton mass ma- 
trices of the model: 

Mu = hKu + fvu Md = hKd + fvd 

= hKu - ^fvu Ml = hKd - ^fVd 

M^f = fvn. (5) 

Here is the Dirac neutrino matrix and is the Majorana mass matrix. 
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Before proceeding, we should specify the origin of CP violation in the 
model. We shall assume that it is spontaneous or soft, that will keep the 
number of parameters at a minimum. The Higgs sector described above 
already has enough structure to generate realistic CP violation either softly 
or spontaneously.^^ The Yukawa coupling matrices h and / in this case are 
real and symmetric. Although there will be three different phases in the 
vev's (one common phase for /t„ and Kd and one each for f„ and Vd)^ only 
two combinations enter into the mass matrices, as the overall phase can be 
removed from each sector. We shall bring these two phases into Vy, and Vd 
and hence forth denote them by fuC*" and f^e*^. 

To see the predictive power of the model as regards the neutrino spectrum, 
note that we can choose a basis where one of the coupling matrices, say /i, 
is real and diagonal. Then there are 13 parameters in all, not counting the 
superheavy scale Vr. 3 diagonal elements of the matrix hKu-, 6 elements of 
fvui 2 ratios of vev's ri = Kd/ Hu and r2 = Vd/vu^ and the two phases a and 13. 
These 13 parameters are related to the 13 observables in the charged fermion 
sector, viz., 9 fermion masses, 3 quark mixing angles and one CP violating 
phase. The light neutrino mass matrix will then be completely specified in 
terms of other physical observables and the overall scale Vr. That would lead 
to 8 predictions in the lepton sector: 3 leptonic mixing angles, 2 neutrino 
mass ratios and 3 leptonic CP violating phases. 

The relations of eq. (5) hold at the intermediate scale M/ where quark- 
lepton symmetry and left-right symmetry are intact. There are calculable 
renormalization corrections to these relations below M/. The quark and 
charged lepton masses as well as the CKM matrix elements run between Mj 
and low energies. The neutrino masses and mixing angles, however, do not 
run below M/, since the right-handed neutrinos have masses of order M/ and 
decouple below that scale. The predictions in the neutrino sector should then 
be arrived at by first extrapolating the charged fermion observables to M/. 

We shall present results for the non-Susy 5*0(10) model with the G224 
intermediate symmetry. We fix the intermediate scale at M/ = 10^^ GeV and 
use the one-loop standard model renormalization group equations to track 



6 



the running of the gauge couphngs between Mz and M/. For Susy-S'O(lO), 
the results are similar, we shall postpone details to a forthcoming longer 

14 

paper. 

To compute the renormalization factors, we choose as low energy inputs 
the gauge couplings at Mz to be 

ai{Mz) = 0.01688 ; a2{Mz) = 0.03322 ; a^{Mz) = 0.12 . (6) 

For the light quark (running) masses, we choose values listed in Ref. (15): 



The top-quark mass will be allowed to vary between 100 and 200 GeV. 

Between 1 GeV and Mz, we use two-loop QCD renormahzation group 
equations for the running of the quark masses and the SU (3)c gauge coupling, 
treating particle thresholds as step functions. From Mz to M/, the running 
factors are computed semi-analytically both for the fermion masses and for 
the CKM angles by using the one-loop renormalization group equations for 
the Yukawa couplings and keeping the heavy top-quark contribution.^^ The 
running factors, defined as 77^ = mi{Mi) / mi{mi) [r]i — mi{Mi)/mi{l GeV) 
for light quarks {u,d,s)\ are r]{u,c,t) = (0.227,0.241,0.452), r]{d,s,h) = 
(0.232,0.232,0.295), r]{e,fx,T) = 0.912 for the case of rut = 130 GeV and 
V{u,c,t) = (0.237,0.252,0.482), r]{d,s,b) = (0.242,0.242,0.302), r]{e,fx,T) = 
0.952 for nit — 150 GeV. The (common) running factors for the CKM angles 
(we follow the parameterization advocated by the Particle Data Group) 5'23 
and Si3 are 1.054 for rrit = 130 GeV and 1.077 for rut = 150 GeV. The 
Cabibbo angle 5'i2 and the KM phase Skm are essentially unaltered. 

Let us first analyze the mass matrices of eq. (5) in the hmit of CP 
conservation. We shall treat spontaneous CP violation arising through the 
phases of the vev's v„e'" and vae^^ as small perturbations. This procedure 
will be justified a posteriori. In fact, we find that realistic fermion masses, 
in particular the first family masses, require these phases to be small. 



m„(l GeV) = 5.1 ± 1.5 MeV 
m,(l GeV) = 175 ± 55 MeV 

mh{mh) 



4.25 ±0.1 GeV 



mail GeV) = 8.9 ± 2.6MeV 
mdmc) = 1.27 ±0.05 GeV 



(7) 
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We can rewrite the mass matrices Mj, and of eq. (5) in terms 
of the quark mass matrices and three ratios of vev's - ri = na/ K'u-i f2 — 
Vd/vu , R = Vu/vr: 

M, = ±^M. - !1±^M. . 

r2 -ri r2 - ri 



D ^ — ' - M. —Ma , 

= l^j_M„-l^— M, . (8) 
Rri-r2 Rri-r2 

It is convenient to go to a basis where is diagonal. In that basis, is 
given by Af^ = yj^f^^fona^yT^ -^j-^gj-g j^diagonai _ (jiagonal(mrf, m^, m^) and 

V is the CKM matrix. One sees that Mi of eq. (8) contains only physical 
observables from the quark sector and two parameters ri and r2. In the 
CP-conserving limit then, the three eigen-vahics of Mi will lead to one mass 
prediction for the charged fermions. To see this prediction. Mi needs to 
be diagonalized. Note first that by taking the Trace of Mi of eq. (8), one 
obtains a relation for ri in terms of r2 and the charged fermion masses. This 
is approximately ri ~ (m^- + 3mi,) /Arrit (as long as r2 is larger than mh/mt). 
Since \mi,\ ~ |m^| at the intermediate scale to within 30% or so, depending 
on the relative sign of mh and m^, ri will be close to either niiy/mt or to 
{mi,/2mt). Note also that if r2 ^ ri, Mi becomes independent of r2, while 
Mjf retains some dependence: 

Mi^AnM^-SMa , M^^M^--Ma . (9) 

r2 

This means that the parameter r2 will only be loosely constrained from the 
charged fermion sector. 

We do the fitting as follows. For a fixed value of r2, we determine ri from 
the Tr(M/) using the input values of the masses and the renormalization 
factors discussed above. Mi is then diagonalized numerically. There will be 
two mass relations among charged fermions. Since the charged lepton masses 
are precisely known at low energies, we invert these relations to predict the d- 
quark and s-quark masses. The s-quark mass is sensitive to the muon mass, 
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the (i-mass is related to the electron mass. This procedure is repeated for 
other values of r2. For each choice, the light neutrino masses and the leptonic 
CKM matrix elements are then computed using the see-saw formula. 

We find that there are essentially three different solutions. A two-fold 
ambiguity arises from the unknown relative sign of rrib and at M/. We 
have found acceptable solutions for both signs. Our numerical fit shows 
that the loosely constrained parameter r2 cannot be smaller than 0.1 or so, 
otherwise the (i-quark mass comes out too small. Now, the light neutrino 
spectrum is sensitive to r2 only when r2 ~ Arris /nic ~ ±0.4, since the two 
terms in M^f of eq. (9) become comparable (for the second family) then. For 
larger values of r2, the first term in M^f dominates and the light neutrino 
spectrum becomes independent of r2- Two qualitatively different solutions 
are obtained depending on whether r2 is near ±0.4 or not. 

Numerical results for the three different cases are presented below. The 
input values of the CKM mixing angles are chosen for all cases to be 5*12 = 
—0.22, 5*23 = 0.052, 5*13 = 6.24 x lO"'^. Since 5km has been set to zero for 
now, we have allowed for the mixing angles to have either sign. Not all signs 
result in acceptable quark masses though. Similarly, the fermion masses can 
have either sign, but these are also restricted. The most stringent constraint 
comes from the rf-quark mass, which has a tendency to come out too small. 
Acceptable solutions are obtained when ^^23; ^13 are in the first quadrant 
and di2 in the fourth quadrant. We use the precisely known charged lepton 
masses at low energies and the running factors discussed above to arrive at 
the values of the masses and mixing angles at M/. 

Solution 1: 



Input : m„(l GeV) = 4 MeV, m^{mc) = 1.22 GeV, rrit = 150 GeV 
mbirrib) = -4.35 GeV, n = -1/51.7, rs = 2.0 
Output : md(l GeV) = 6.4 MeV ms{lGeV) = 164 MeV 




9 



lepton 



KM 



0.9522 
-0.2991 
-0.0615 



0.3051 
0.9400 
0.1522 



0.0123 
-0.1637 
0.9864 



(10) 



Solution 2: 

Input : m„(l GeV) 

rribimb) 
Output : md{l GeV) 



V- lepton 
KM 



Solution 3: 

Input :m„(l GeV) 

mbirrib) 
Output : m<i(l GeV) 



V 



lepton 



KM 



4 MeV, mc{mc) = 1.22 GeV, rut = 150 GeV 
-4.35 GeV, n = -1/51.4, r2 = 0.2 

5.6 MeV, m,(l GeV) = 175 MeK 
R (6.3 X 10"^ 2.1, -2.7 x lO"^) GeV 

0.9969 0.0506 -0.0607\ 
-0.0585 0.9890 -0.1359 . (11) 
0.0532 0.1390 0.9889 / 



3.5 MeV, mc{mc) = 1.27 GeV, rrit = 130 GeV 
-4.35 GeV, n = -1/101.8, r2 = -0.5 
-5.24 MeV, m,(l GeV) = -173 MeV 
R (4.1 X 10"^ 1.1, 6.2 X 10^) GeV 

/ 0.9954 -0.0794 -0.0521 \ 

0.0704 0.9853 -0.1559 . (12) 
V 0.0637 0.1515 0.9864 / 



Solution 1 corresponds to choosing ri ~ mi,/ nit . All the charged lepton 
masses are negative in this case. Since r2 is large, the Dirac neutrino matrix is 
essentially M„, which is diagonal; so is the Majorana matrix. All the leptonic 
mixing angles arise from the charged lepton sector. Note that the predictions 
for rrid and are within the range quoted in eq. (7). The mixing angle 
sin^,^^_,^^ relevant for solar neutrinos is 0.30, close to the Cabibbo angle. Such 
a value may already be excluded by a combination of all solar neutrino data 
taken at the 90% CL (but not at the 95% CL).^ Actually, within the model, 
there is a more stringent constraint. Note that the — u-,- mixing angle is 
large, it is approximately S\Vcb\ — 0.16. For that large a mixing, constraints 
from 1/^ — Vt- oscillation experiments imply^^ that |m^^ — ml \ < 4 eV"^. 
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Solution 1 also has m^^/m^^ ^ 2.5 x 10^, requiring that m^^ < 0.8 x 10~^ eV . 
This is a factor of 2 too small for z/g — MSW oscillation for the solar puzzle 
(at the 90% CL), but perhaps is not excluded completely, once astrophysical 
uncertainties are folded in. If Vr mass is around 2x 10^'^ eV , Ve — ^r oscillation 
may be relevant, that mixing angle is ^ ^\Vtd\ — 6%. It would require the 
parameter R = Vu/vr ~ lO^"*^^ or vr ~ 10^^ GeV for f „ ~ 1 GeV. Such a 
scenario fits well within Susy-S'O(lO), but not in the non-Susy 6*224 chain. 
Note that Ur mass is negative, a transformation Ur iVr will make it positive. 

Solution 2 differs from 1 in that r2 is smaller, r2 = 0.2. The 1 — 2 mixing 
in the neutrino sector is large in this case, so it can cancel the Cabibbo like 
mixing arising from the charged lepton sector. As we vary r2 from around 0.2 
to 0.6, this cancellation becomes stronger, the — mixing angle becoming 
zero for a critical value of r2- For larger r2, the solution will approach Solution 
1. The — Vr mixing angle is still near 3|Vcfe|, so as before, m^^ < 2 eV . 
From the v^-|v^;^ mass ratio, which is 1.3 x 10^ in this case, we see that 
'niy^ < 1.6 X 10~^ eV. This is just within the allowed range'^ (at 95% CL) for 
small angle non-adiabatic — MSW oscillation, with a predicted count 
rate of about 50 SNU for the Gallium experiment. Note that there is a lower 
limit of about 1 eV for the Vr niass in this case. Forthcoming experiments 
should then be able to observe — Vr oscillations. A mass in the (1 to 2) 
eV range can also be cosmologically significant, it can be at least part of the 
hot dark matter. In Susy 5*0(10), Vf. — Vt oscillation (the relevant mixing is 
about 3|Vt(i| — 5%), could account for the solar neutrino puzzle. 

Solution 3 corresponds to choosing r\ ~ {mi,/2mt). All charged lepton 
masses are positive in this case. The sign of r2 has been chosen to get 
small sin^e-/i- (For other values of r2, the results are similar to Solution 1.) 
However, the mass ratio z^r/^'^t is ~ 6 x 10^, and sin^^^^ ~ "SVcb so z/g — 
oscillation cannot be responsible for solar MSW. As in other cases, u^. — 
MSW oscillation with a 6% mixing is a viable possibility. 

Observe that none of the solution generates Vf. — z/^ mixing large enough 
for the vacuum oscillation for solar neutrinos. Similarly, the puzzle with 
atmospheric neutrinos cannot be explained in this minimal scheme in terms 
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of Vp^—Vr oscillation, the relevant mixing is not large enough. (For an 5*0(10)- 
based explanation of this phenomenon, see Ref. 18). 

Let us now re-instate the CP-violating phases a and /5 in the vev's per- 
turbatively. Small values of the phases are sufficient to account for realistic 
CP violation in the quark sector. We shall present details for the case of So- 
lution 2 only, others are similar. We also tried to ffi all the charged fermion 
masses and mixing angles for large phases, but found no consistent solution. 

First we make a basis transformation to go from the basis where M„ is 
diagonal to one where the matrix Hku is diagonal. It is easier to introduce 
phases in that basis. For a = 3.5*^, (3 = 4.5°, the CP-violating parameter J 
for the quark system^^ is J ~ 1 x 10~^, which is sufficient to accommodate e in 
the neutral K system. The leptonic CP violating phases are correspondingly 
small, for eg., the analog of J is J; ^ 7 x 10^^. These small phases modify the 
ffist family masses slightly, but the effect is less than 10%. Our predictions 
for the neutrino mixing angles are essentially unaltered. 

In summary, we have presented a class of minimal 5*0(10) grand unified 
models where the light neutrino masses and mixing angles are predictable 
in terms of observables in the charged fermion sector. Our approach here 
has been orthogonal to some other recent attempts based on grand 
unification, we have kept the Higgs sector as simple as possible and followed 
its consequences. We have found three different types of solutions for the 
neutrino spectrum. In Solution 1, the I'e—'^^i mixing angle is near the Cabibbo 
angle, while Solutions 2 and 3 have it much smaller. In all cases, z/g ~ '^t 
mixing angle is predicted to be near 3\Vtd\ — 0.05 and — mixing angle 
is ~ 3|Vcfe| — 0.15 with the mass ratio m^^/niy^ > 10^. If the solar neutrino 
puzzle is due to small angle non-adiabatic MSW, as in Solution 2, t'^ — t^r 
oscillation should be observable in the forthcoming experiments. 
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